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1. INTRODUCTION 

Laser ranging systems have developed to the point where it may now be 
possible to measure the distance from an orbiting satellite to a point on the 
earth with uncertainties on the order of only a few centimeters. The precision 
of satellite ranging systems will be limited in part by atmospheric refraction 
and scattering. In this paper the effects of atmospheric turbulence on the 
accuracy of single color and multicolor ranging systems is discussed. The 
statistical characteristics of the random path length fluctuations induced 
by turbulence are examined. Correlation and structure functions are derived * 
using several proposed models for the variation of C^ with altitude. For single 
color systems it is shown that the random path length fluctuations can limit 
the accuracy of a range measurement to a few centimeters. Two color systems 
can partially correct for the random path fluctuations so that in most cases 
their accuracy is limited to a few millimeters. However, at low elevation angles 
and over long horizontal paths two color systems can also have errors approaching 


a few centimeters. 
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2. RANDOM PATH LENGTH FLUCTUATIONS 
2 . 1 Theoretical Background 

Refractive Index fluctuations caused by atmospheric turbulence will 

induce random variations in the optical path length. The magnitude of these 

path length variations can be estimated from the theoretical analyses of 

optical phase fluctuations published by Tatarski [1] and other workers. The 

instantaneous path diviation AL can be predicted using geometric optics 

• AL = /dr n^r) . (1) 

. C 

nj_ is the fluctuating part of the refractive index about the mean value and 
C is the ray path. If random variations in the ray path due to n^ are neglected,* 
C will be determined by the average refractivity. The average path deviation 
is zero because the average value of n^ is zero 

<AL> = /dr <n 1 (r)> = 0. (2) 

C 

The covariance between the deviations along two different paths can be 
written in terms of the refractive index covariance 

b at = ’/- d Ii/ d ^2 <niCri)ni(r 2 ) > . (3) 

Ci C ? 

We will assume that the refractive index is isotropic and its covariance is 
given as the product of a function depending on the magnitude of the difference 
coordinate and another function depending on the average coordinate 

rj + r ? 

< n 1 Cri)n 1 (r 2 )>=B^C 2 (4) 

B° measures the correlation between the fluctuations at points r x and r 2 while 
measures the strength (Variance) of the fluctuations. When the paths are 
identical (C^ = C 2 ) , C3) gives the mean-square deviation. 

* It can be shown that variations in C due to n will increase the optical path 
length by at most a few nm. 1 
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Two different path, geometries will be considered in this paper (Figure 1) . 
The angular- path separation would be encountered with a stationary ranging site 
and a moving target such as a satellite. The parallel path separation would be 
encountered when the ranging site and target are moving in the same direction 
at identical speeds or when both site and target are stationary and there is a- 
uniform wind blowing transverse to the path. To simplify the problem we will 
assume that the propagation paths are straight lines. Refractivity gradients 
in the atmosphere will cause the actual propagation paths to deviate from a 
straight line by as much as a few meters, particularly along the satellite path 
at low elevation angles. However, the deviation will only be significant at 
the higher altitudes where is small. Hence, the error will also be small. 
For horizontal paths, the curvature is negligible. 

The parallel path covariance is calculated by substituting (4) in (3) and 
noting the geometry in Figure 1 

L L j 

B AL (d) = 4 / dr x /dr 2 B° [((rj - r 2 ) 2 + d*)*] . (5) 

o o 

The constant 4 appears in equation (5) because we are considering a complete 
"roundtrip" between the ranging site and target. Equation (5) can be rewritten 
in a more convenient form by making the change of variables (rj + r 2 )/2 = £ and 
Oi - r 2 ) = p 

L L-p/2 ' , 

B.t (d) = 8 /dp /d? B^CO B° Up 2 + d 2 ) 2 ]. (6) 

AL o p/2 n n 

The correlation length of the refractive index fluctuations is on the order of 
the outer scale of turbulence, L . If the path length L is much larger than 
L q , the upper limit on the p integration can be set equal : to infinity and the 
5 integration can be extended from 0 to L 

CO 

B AL (d) = 8 /d? B^CO /dp B° [(p 2 + d 2 fy. 
o o 


( 7 ) 
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TARGETS 


PARALLEL PATH GEOMETRY 


Figure 1. 


Geometry of the laser ranging site and targets. 



. 5 


The error introduced by this approximation is on the order of . 

Lt 


can be conveniently described in terms of its spatial spectrum 


B°(r) = 4r /dK K 2 $° (K) S1 ^ r:) . 
n J n Kr 

o 


( 8 ) 


The modified von Karman spectrum is commonly used for $ 

-K 2 /K 2 
ae m 

$°(K) = * r-_ 

n • (1 + k 2 l 2 ) 11/6 


n 


K = 5.92/1 (9) 

m o 

L o rCii/6) 
a r(i/3) 


where a is a constant chosen so that' B°(0) = 1 and l is the inner scale of 

n v o 

turbulence. The von Karman spectrum appears to be a.dequate under conditions 
where the low frequency behavior of the spectrum is unimporant. However, the 
random path variations are predominantly influenced by- refractive index 
perturbation's in the spectral input range (K<2ir/L o ) . Therefore, our results 
will be developed in terms of an arbitrary so that the effects of different 
spectra can be studied. Substituting (8) in (7) and carring out the p integration 
gives 


L 00 

B AL (d) = 16m 2 /d£ B^(0 /dK K$°(K) J (Kd) . (10) 

o.o ' 

Another statistic of interest is the structure function for the path 
deviation. It is defined as the mean-square path difference and can be 
written in terms of B,. 


D AL (d) = 2[B AL (0) - B AL (d)] = 32 it 2 /d£ B^(§) /dK K$°(K) [1 - J„(Kd)l 


n 


(ID 
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When the path separation is small, the structure function is not sensitive to 
the low frequency behavior of 4^. 

Similar expressions can be derived for angular path separations. Noting 
the geometry in Figure 1 and substituting (4) in (3) we obtain 

B {-a) = 4 /dr /dr B V cos C<l>/2) 3 B°[((r -r ) 2 cos 2 (4>/2) + (r +r ) 2 sin 2 (<|>/2))' S ] - 

AL o 1 o 2 n Z n 1 2 12 

( 12 ) 

Equation (12-) was derived from (4) by assuming that is a function of the 
average of the projections of rj and r 2 onto the bisector of the angle <|>. By 
making the change of variables 


apices M/2) = 5 
(r 1 -r 2 )cos(<f>/2) = p 


(13) 


and extending the p integration to infinity we have 

L “ j 

B AT (40 8 /dg B^C?) /dp B°[(p 2 + 4C 2 tan2( (f) /2))1]. (14) 

cos z (*/2) o n o n 

It is most convenient to write the covariance in terms of the refractivity 
spectrum. -Substituting (8) in (14) and integrating over p gives 

B AL<« " cos‘l!/2) ^ B n ( « J 0 [2K«tan(t/2U. (15) 

o o 

The structure function for the path deviation is 

L ■ J [2K€tan($/2)] 

• D C«-32,r* /dKK*°00 /« BJJCO [1 - ° C0S 2 ( ^ --1. (16) 

o’ o 
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The path deviation statistics for both parallel and angular path 
separations are summarized in Table I. These statistics were derived assuming 
that the refractive index correlation length L q was much smaller than the 
total path length L. In this case the path deviation given in equation (1) is 
the sum of many statistically independent values of the refractive index. The 
central limit theorem can be used to argue that AL is Gaussian distributed. 
Thus, the correlation functions in Table I provide a complete statistical 
description of the Gaussian random variable AL for the path geometries 
considered. In the following sections these statistics ' will be evaluated for 
typical atmospheric parameters . 
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Covariance 


Structure- 

Function 


Covariance 


Structure 

Function 


TABLE I 

' i 

Path Deviation Statistics 


Paralle.l. Path 


L 

B AL Cd) = 1,6tt 2 Jd? B V (£) / dKK$ 0 (K) J (Kd) 
o • o 11 

L " 

D AL (d) = 32ir 2 /d? B^ (O /dKK$°CKHl - J Q (Kd).] 
o o 


Angular Path 


16tt 2 


/dKK*°(K) / d? B^C?) J [2K5 tan (*/2-) 3 


T C*fO - * v " J ' " " rt * 

“ J-* 9 /i' n n. o 

cos^(<j>/2) o o 


t> AL C<j>) = 32tt 2 /dKK^CK) /de B^C5) [1 - — 


J (2K? tan (<j>/2)) 


n' n 

o o 


cos 2 C4>/2) 


-] 
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2.2 RMS Path Deviation and Correlation Function 
The mean square path deviation can be calculated using the covariance 
functions by setting d or <j> equal to zero in (10) or, (15) 

L °° 

<AL 2 > = 16ir 2 / d£ B V (?) /dfC K4°(K). (17) 

o o 

The refractive index variance, B v , is related to the turbulence structure 

n 

parameter [ 2] 

B V (r) = 0.033C 2 ( r ) L 2/3 = 0.523 C 2 (r)L 2/3 . (18) 

n n ' r(ll/6) o no 

It is convenient to write the path deviation in terms of the structure parameter 

at the' ranging site and an effective path length, L , defined as 

0 

L e - ld( ‘S® ' C19) 

n o 

The mean-square path deviation for the von Karman spectrum (eq. (9)) is 

' t 

<AL 2 > = 0.033 C 2 (o) L 5/3 L ~ =~ 3. 127 C 2 (o)L 5/3 L 

n 5 o -e n o e (20) 

The rms path deviation is plotted in Figures 2 and 3 as a function of C 2 for 

several values of L and L . 

e o 

For horizontal paths the structure parameter will be essentially constant 

so that L is equal to the average path length L. Measurements by Bufton [3,4] 

0 

have shown that turbulence exhibits a layered- structure as a function of altitude. 
Except for a sharp peak near the tropopause, the structure parameter also exhibits 
an overall decay with increasing altitude (Figures 4 and 5) . This decay roughly 
follows an exponential model [4-6] 

C 2 (h) = C 2 (o)e" h/h s • (21) 



L 0 =50m 10m Im 


IO" IS IO' 14 IO' 13 IO' 12 

C 2 (m“ 2/3 ) 


Figure 2. Rms path deviation given by eq. (20) for the von Karman spectrum. is the 
refractive index structure constant, L 0 is the outer scale of turbulence and 
L e is the effective path length in turbulence defined by eq. (19). 




L 0 =500m 200m 100m 


C„ (m l 


Figure 3. Rras path deviation given by eq. (20) for the von Karman spectrum. C 2 is the 
refractive index structure constant, L Q is the outer scale of turbulBncfe and 
L e is the effective path length in turbulence defined by eq. (19). 
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from Bufton's report, ref. [3]. 



Figure 5. 'Hufangel's updated model for (^(h): , sunny day; 

; clear night; dawn-dusk minimum; . .., 

disturbed layers. Local ground is assumed level and 
not on a mountain top. Figure is from Brookner's 
paper, ref. [5]. 
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where h g is the atmospheric scale height and h is the altitude. h s is typically 
on the order of 2 to 5 Km. Then for satellite ranging the effective path 
length is 


L = h csci|> (22) 

e s 

where is' the satellite elevation angle. L q also varies with altitude. Values 

frequently used for L are about 100 meters or 1/3 to 1/5 the height above ground, 

whichever is less. For satellite ranging 100 meters .should be used for the 

value of L . 

o 

The path deviations are predominantly influenced by the input, range of 
the refractivity spectrum (K^tt/L^) . To investigate the sensitivity of the path 
deviations to the form of the mean-square deviation was calculated for the 


following spectrum 


*°(K) 


-k 2 /k 2 

e m 

(1+K 2 L 2 ) V 

o 


(23) 


<AL 2 > 


r(i/3)r(yt-i) 

r(5/6)r(p-3/2) 


[0 . 033C 2 (0) 
n 


48v 2 

5~~ 


l 5/3 

o 



•= 7.42 


r(u-i) 

r(y-3/2) 






(24) 


The ratio of the rms deviation for arbitrary y to the deviation for y = 11/6 is 
plotted in Figure 6. As y increases., more of the energy is concentrated in the 
low spatial frequencies for a fixed value of the refractive index variance. This 
effectively lengthens the tail of the refractive index covariance function which 
increases the rms deviation. 

The parallel j>ath correlation function can easily be calculated by substituting 

(23) into (10) and integrating over K. If inner scale effects are neglected (i.e. 

let K =co in (23)) the correlation function can be evaluated in closed form 
m 


R AL< d > * 


B AL Cd)/B AI.™ = 


2 2 -M 

r(u-i) 




(<L> 

•1 '■L J 
o 


(25) 


where K , is a modified Bessel function. 
y-1 






The correlation functions are plotted in Figure 7 for several values of y. 
When y increases, the correlation tails increase as expected. The deviations 
on the two paths are essentially uncorrelated when the' separation is on the 


order of L . The effects of the inner scale of turbulence are illustrated in 
o 

Figure 8 where the correlation function is plotted for y = 11/6 for several 

values of K L . K L is typically greater than 10 3 ' in the atmosphere so that 
mo, mo r 

inner scale effects can be neglected. 

The angular path covariance function is more difficult to evaluate because 
it depends on the detailed behavior of the atmospheric structure parameter. For 
satellite ranging we will use the exponential model for C 2 given in (21) . 
Brookner [5] has proposed a more complicated form of this model which includes 
a component to describe the tropopause contribution. Unfortunately the use 
of his model leads to considerable mathematical difficulties in calculating 
the angular correlation function. For this reason we will use the simpler 
model given in (21) which does give a fairly accurate description of C 2 [4,5]. 

In the following section Brookner's model is used to study the path structure 
function and is compared with the model in (21) . 

Using (21) the ' angular .covariance function is given by 
C3(0)h 


B at (<j>) =' /dK K$°(K) [1 + K 2 p 2 taii 2 (cj>/-2) ]-~ 

AL . cos 2 (<j>/2) o n ' ■ S 


(26) 


where 


p s = 2h s c sc(^+<])/2), 


(27) 


is the initial elevation angle of the satellite (See Figure 1) . 

If the spectra given in equation (23) are used, the angular correlation 
can be written in terms of a hypergeometric function [-7] 


1 fu-11 

R AL^ = B al ( 0)" = cos 2 (4/2) (y-1/2) 


2 Fi(h, 1, U+k, I" tan 2 (<j>/2)) 


(28) 



Figure 7, 


Parallel path correlation function for the spectral 
model in eq. (23) . L 0 is the outer scale of turbulence 
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0.00 1.0 2.0 3.0 4.0 5.0 


d/L 0 


Figure 8. 


■ Parallel path correlation function for the von Karman 
spectrum ( y = 11/6). 
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The path deviations are essentially uncorrelated when 

L sint/j L sini|) 

<j>=2tan _1 (-|j 0 “ = <P • (29) 

s s s 

The correlation angle 4> s is at most only a few degrees. ' When <J> is on the 
order of or less than <f> s , the cosine factor in (28) can be set equal to one 
and the tangent approximated by the first term in its Taylor series expansion. 

In this case the angular correlation function becomes 

R al (« =^qf2?i(h, I, v+k, 1-<P 2 U 2 ) (30) 

L sinili 
^ _ o 


The correlation functions given by (30) are plotted in Figure 9 for y=3/2, 

11/6, and 5/2. When y increases the correlation tails also increases,. 

It should be pointed out that the path fluctuations on the satellite- 
earth path are nonstationary. Because the mean-square path deviations depend 
on the satellite elevation angle (see (20) and (22)), the rms -path .deviations 
will change as the satellite passes over the ranging site. The correlation 
function also depends on the satellite position since and <J> are functions 
of t . The path deviations decorrelate within a few degrees of angular 
separation. In the region of small <j>, p g will be essentially constant so that 
the effects of nonstationarity on the correlation function are negligble. 

For large <j> (>10°) , the nonstationary effects must be considered by using 
equations (26) and (27) or (28) to calculate the correlation. 
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2. 3 Structure Functions 

The structure functions can be evaluated directly from the covariance 
by using the relation in equation (11). In this section we will compare the 
structure functions calculated using the von Karman spectrum with functions 
calculated using a new spectral model proposed by Greenwood and Tarazano [8]. 

The new model was derived empirically from atmospheric microtemperature 
measurements by varying parameters until a best fit to the tempefature 
structure function data was obtained. Greenwood and Tarazano' s (G-T) 
proposed spectrum is given by 

. 3>°(K)<^(K 2 l 2+KL o ) _11/6 . (32) 

It is plotted in Figure 10 along with the von Karman spectrum. 

It is most convenient to calculate the normalized structure function 
because of its relationship to the correlation functions 

o AL ( d ) • 

D -(») ' 1 " R AL^* 

AL 

When d is small compared to the correlation length, the structure function 
is relatively insensitive to the low frequency behavior of the refractivity 
spectrum. In this region - the structure function yields detailed information 
about the behavior of the main lobe of the correlation function. This information 
can be used to describe angle- of- arrival fluctuations and path variations across 
the receiving telescope aperture. In multicolored ranging systems, each 
frequency will traverse a slightly different path to and from the target 
because of dispersion in the atmosphere. The path fluctuations at each 
frequency will also be slightly different. The structure functions can be 
used to estimate these differences so that their effects on ranging accuracy 
of multicolored systems can be analyzed. 




■posed by Greenwood and Tarazano [ 
Karman spectral model, eq. (9]. 
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Using the von Karman spectrum* the parallel-path structure function is 
given by 


D. T (d) ,1/6 , , 

AL k _ , _2 fdj 5/6 ,d_, 

F(5/6) h/6\ } ' 


(34) 


The asymptotic behavior for d small with respect to L q is 


P Ai/ d<< V 

D AL W 


rq/6) 


5/3 (f-) 5/3 = 1.864(f-) 5/3 . 

2 5/3 r ('ll/e) L o L o 


(35) 


When d » L q the asympotote is of course 1. The intersection of D^(d « L q ) 


and D (d « L ) is at 
AL o 


d = 0.688L . 

o 


(36) 


The structure function is plotted versus d/L Q in Figure 11. 

When the G-T model spectrum is used, the structure function cannot be 
expressed as simply as before. Following Greenwood and Tarazano's [8] 
approach we write D in terms of confluent hypergeometric functions 

Zl Li 


D 


AL 


AL 


(d) 

(“T 


in - j 

iT/de U(l/6, -2/3, -i~cos0) 


U(l/6, -2/3,0) 


(37) 


The derivation of (37) is detailed in Appendix A where the first term of the 
power series expansion in j— is also evaluated. The small d asymptote is given by 


D AT (d « L ) 


AL o J - 0. 2093 (y— ) 


d ,5/3 


D AL« 


and the intersection of D (d « L ) and D (d » L ) is at 

Zl J-J O- Z\ Lr O 


(38) 


d = 2.56L . 

o 


( 39 ) 




d/L 0 


Figure 11. Parallel path structure functions, (a) von Karman spectrum, 
eq. (34). (b) Greenwood and Tarazano spectrum, eq. (37). 
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Equations (38) and (39) are considerably different from the corresponding 
results for the von Karraan spectrum ( (35) and (36) ) . This is illustrated in 
Figure 11 where the structure functions for both the von Karman and the G-T 
spectra are plotted. The structure function increases much more slowly with 
the G-T model. This indicates that the corresponding correlation function has 
a much longer tail. The effect results because there is more energy in input 
range of the G-T model (see Figure 10). 

These results are interesting because of their implications regarding 
the mean-square path deviations. The mean square deviation calculated using 
the G-T model is given by 


Ai2> = . J* C^/3) r q 033C^(0) L 1 

G-T T(5/6) L U * U ^V J 5 o e J 


(40) 


= 26.31 C 2 (0) L 5/3 L 
n o e 


The rms deviation calculated for the G-T model is approximately 3 times greater 
than the corresponding results for the von Karman spectrum. Because the 'G-T 
model was derived emperically from actual atmospheric measurements we believe 
that (40) represents a more realistic value for the mean-square path deviations 
than the results obtained from the von Karman spectrum. The rms deviations for 
the G-T model can be obtained from Figures 2 and 3 simply by multiplying the 
values by 


2A (1/3) ,3s 
t T(5/6) J ~ 


2.9 


( 41 ) 



25 


The structure functions for the angular path depend on both the refractivity 
spectrum and the behavior of If is described by the exponential model 

given in (21) , the structure function can easily be calculated from (30) with 

V = 11/6 


d alW 


~ 1 - 


s 

821 


F (h, 


i, 7/3, 1 - <r/<^). 


(42) 


The small angle asymptote is 

^ H . nx_/6jr(4/a . 2 . 804( |_ ) 

V 0 rft) V 

and the intersection of D^ L ((j> <<: and D AL (<i> >:> $ 3 ) at 

<f> = 0 . 539 (J> . 

T 


5/3 


(43) 


(44) 


Equation (42) is plotted in Figure (12) . 

When the G-T spectrum is used, the angular structure function can be 
derived from the parallel path structure function given by (37) - The calculations 
are detailed in Appendix B and plotted in Figure 12 along with the corresponding 
results for the von Karman spectrum. The small 4> asymptote is 


D. r (<j> <J> ) A C/Z 

■ AL n f - V S - 0.315 (^-) 5/3 
D AL W 4 s 


(45) 


and the intersection of D^(<i> << and >> <J> S ) is at 

4 = 2.00 4 . (46) 

The angular structure functions also increase more slowly when the G-T model 
spectrum is used. These results imply that the path deviations do not decorrelate 
as rapidly as might be expected with the von Karman model . 

Brookner [5] has proposed a more complicated version of the exponential 
model for the variations in with altitude 



u> ,v a/(<£) 
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Figure 12. Angular path structure functions. (V) von Karman spectrum, 
eq. (42). (b) Greenwood and Tarazano spectrum, eq. (45). 
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C £< h > = C L h " b e " h/hs + C nT 5 < h -V (4 

h is the atmospheric scale height which in general is different from the h 

s 

employed in (21) . The parameter b is chosen depending on the conditions, and 
C^Q specified so that C 2 (h) is in m when h is in meters. The second 
term in (47) is the tropoause contribution. h T is the tropopause height and 

.j. ‘ 

C nT is the integral of C 2 (h) under its peak at the tropopause. When b and 

C* T are zero, (47) reduces to the model given in (21). Brookner suggests a 

value of 320 m. for h and 12 km. for h . The model given by (21) requires h 

Si s 

to be on the order of 3. 5 km. 

When Brookner' s model is used, the covariance and structure function 

separate into two terms. To illustrate, the angular correlation function 

is derived in Appendix C for the von Karman spectrum using Brookner' s model 

for C 2 (h) 
n 


R 4L C « 


n 


+ + 
C C n t( W 


B ' F ft-b/2, l-b/2, 7/3-b, l-* 2 /*?] 
° 2 1 s 


~ n J _- g (<j>/<j> ) 5/6 JC , (*/<}> ) , b<l 

KKt> 1 t 5/6 


(48) 


where 



L 

♦t = 4 sin * 


(49) 


+ 4 - 

and 8^ and B are constants defined in Appendix C. and are the integrated 
contributions due to the lower atmosphere and tropopause respectively. In this 
case th.e path deviations are characterized by two distinct correlation angles. The 
tropopause correlation angle ^ is generally much smaller than the lower atmosphere 
correlation angle 4> s - Consequently, there is an initial rapid decorrelation as <j> 
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approaches (j)^. The correlation then saturates at a value approximately 
equal to C /(C^C^) until 4> approaches 4> g . Then the correlation decays 
to zero as the deviations due to the lower atmosphere decorrelate. 

This effect is nicely illustrated in Figures 13 and 14 where the 
correlation and structure functions are plotted. We have used the parameter 
values suggested by Brookner in reference [5]. The values were selected by 
Brookner to obtain a best fit of the model to data published by Hufnagel [9]. 
In all cases the initial decorrelation (increase of the structure function) is 
dominated by the tropopause contribution. 

It should be pointed out that even Brookner 's improved model is 
still an approximation at best. Bufton's [3, 4] data which was collected 
during the dawn-dusk minimum periods show a highly layered structure to the 
turbulence. Based upon the results from Brookner* s model, we might expect 
the strongest layers to each have an associated correlation angle so that 
the structure functions would be even more complicated than those illustrated 
in Figures 13 and 14. The initial decorrelation will probably be determined 
by the highest layers near the tropopause. 

The effects of turbulence on laser ranging systems will be 

most apparent on sunny days when the turbulence is strong. The G-T 

spectrum predicts rms deviations approaching a few centimeters (eq. (40) 

and (41) arid Figure 3). In this case the lower atmosphere is the major 

source of the path deviations. Significant decorrelation will occur 

only when the separation approaches the lower atmosphere correlation angle 

a while the initial decorrelation will still be determined by the 
s 

tropopause (Figures 13 and 14) . 



Figure 13. Angular path correlation function calculated using Brookner's [5] 
improved model for Cr (h) given in eq. (47). (a) Sunny day. 

(b) Clear night, (c 1 ? Dawn-dusk minimums. 



Figure 14-a. Angular path structure function calculated using Brookner's [5] 
improved model for C2 (h) given in eq. (47). (a) Sunny day. 

(b) Clear night. (cj Dawn-dusk minimums. 
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Figure 14-b. 


Angular path structure function calculated using Brookner's [5] 
improved model for (h) given in eq. (47). (a) Sunny day. 

(b) Clear night. (cj Dawn- dusk minimums. 
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3. EFFECTS OF PATHS FLUCTUATIONS ON RANGING ACCURACY 

The statistical characteristics of random path length fluctuations were 
derived -in the previous section. In this section the effects of path fluctua- 
tions on the accuracy of single color ^and two color ranging systems are 
examined. 

3. 1 Single Color Ranging Systems 

The optical path length for a roundtrip between the ranging site and 
target is given by 

L = /di^n (r) + Jdr n^r) = <L> + AL (50) 

c c 

<L> is the average path length which is determined by the average refractivity 
(n Q ) along the ray path. AL is the instantaneous random path deviation caused 
by fluctuations in the refractive index (nj) . AL is Gaussian distributed and 
its statistics were derived in Section 2. 

When a range measurement is made, L is the parameter that is actually 
measured. Formulas are available (e.g. ref. [10-12]) which can partially 
correct the measurements for the effects of the average atmospheric refraction. 
Since AL is random, it is not possible to correct the measurement for the 
effects of turbulence. AL represents a fundamental limit in the accuracy of 
single color ranging systems. 

The rms value of AL was calculated in the previous section and plotted •. 
in Figures 2 and 3 for typical atmospheric parameters. Figures 2 and 3 are 
based on equation (20) which was derived for the von Karman spectrum. These 
results probably represent a conserative estimate for AL . The rms 
deviation was also calculated using a new spectral model proposed by Greenwood 
and Tarazano [8]. The G-T model gives a value for approximately 3 times 

larger than the results for the von Karman spectrum. Because the G-T model 



was derived emperically from actual atmospheric measurements of the temperature 
structure function, we believe that it gives a more realistic value for the 
rms path deviations. 

When the structure parameter follows the exponential model in (21) the 

effective path length in turbulence is given by (22) . L q and h g -are on 

the order of 100 m and 3 km respectively. Under these conditions the rms path 

deviations can be up to a few centimeters when the satellite is at low elevation 

angles (-10°) and the turbulence is very strong (C 2 -10 -13 m -2 / 3 ). Under most 

conditions- C 2 will be much .weaker (10" 15 m -2 / 3 ) so that rms deviations will 
n 

be a few millimeters or less. 

It is possible to reduce the effects of the path deviations by averaging 
many different range measurements. The amount of reduction depends on the 
correlation between the path deviations on each measurement. If the path 
deviations are stationary and n measurements are averaged, the residual rms 
deviation is given by 


1 


n 


n 


( 2 

i=l 


n 

l 

j=l 




AL 


rms 


( 51 ) 


where AL is the deviation on a single measurement and p . is the 
rms 

correlation between the path deviations on i^ 1 and measurement. If the 
measurements are statistically independent, the rms deviation is reduced by 
the factor l//n~\ 

The path deviation correlation functions were derived in Section 2 (see 

Figures 7 and 9 ). For satellite ranging the path deviations decorrelate 

L 0 

when the angular separation is on the order of <j> s = -g— sin'f'. The time 

s 

required for the Lagoes satellite to move an angular distance is on the 
order of tens of seconds. Typical laser ranging systems are capable of 
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making many measurements within this time period. When the correlation between 
successive measurements is high, the summations in (51) can be approximated 
as integrals. ,To see this, let the angular velocity of the satellite as viewed • 
from the ranging site be constant and denoted by w. The averaging operation 
can be written as an integration 

1 T 

AL avg = i / dt AL(<j>=i|H-a>t) (52) 

T is the length of the time interval over which the measurements are averaged 
and ^ is the satellite elevation angle on the first measurement at t=o. The 
mean square value of the integral is given by 

• • - x 

<ALj vg > = | /dt (l-t/T) B^ l (<J> = tot) (53) 

Equation (53) can be rewritten in terms of the refractive index spectrum 


<AL 2 >' 
avg 


64ir 2 

cjT 


L tan(mT/2) 

JdK K$°(K) / d? B^(0 /dx(i-~(tari' 1 X) J Q .(2Ka) . 

O 0 0 


(54) 


(54) was derived by substituting (15.) into (53) and making the change of 
variables x = tan (wt/2) . Notice wT is the total angle over which the measure- 
ments were averaged. 

The ratio of Al for the averaged measurements to AL for a single 
rms rms 6 

measurement is plotted in Figure 15 as a function of toT. The results were 
calculated using the angular correlation function given in (28) for p = 3/2, 
11/6 and 5/2. This corresponds to the spectral model given in (23) for 9>° 
and the expontential model given in (21) for B^. There is suprisingly little 
reduction even when the measurements are averaged over 10 correlation angles. 
This is due primarily to the long tails on the angular correlation functions 
(see Figure 9) . The effects of iidnd were not included in the derivation 
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Figure 15. Ratio of the rms path deviation for time-averaged measurements 
to the rms deviation for a single measurement, to is the 
angular velocity of the satellite as' seen from the, ranging 
site^ T is the averaging time and <f> is the correlation angle 
(eq. (30)). The results we calculatld using the spectral- model 
given by eq. (23). 
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of the correlation functions in the previous section. Wind motion transverse 

i 

to the propagation path will tend to decorrelate the path deviations more 
rapidly when the wind velocity is in the direction opposite the satellite motion. 
When the wind velocity is in the same direction as the satellite motion, the 
path deviations tend to remain correlated much longer.. Figure 15 was derived 
assuming no wind and probably represents an average of these two extreme cases. 
Wind motion will be important only when the transverse speed Vj_ is large 
enough to blow a turbulent blob of size L q across the propagation path during 
the integration time, i.e. if T > L o /Vj L . The wind speed in the’ lower atmosphere 
is typically on the order of a few tens of meters per second. Therefore, the 
results in this section apply for integration times from a few seconds to about 
1 minute . 

Equations (52) and (54) can be simplified if we assume that the integration 
time is long compared to the correlation time (i.e. coT»<j> s ) . In this case t/T 
will be small over the important range of integration and can be neglected in 
(52). The upper limit on the T integration can be extended to ir/tu since B 

, Li 

is negligible when T»^ s /w 

ir/w 

<AL 2 > - i fdt B AL (<j.=wt) . (55) 

& o. 

Similarly, if we neglect the dependence of on 4“ (B^ is evaluated for the 
elevation' angle $ + <f>/2) , the term tan^X can be neglected in (54) and the 
upper limit on the X integration can be extended to infinity 


<AL 2 > - 
avg 


32ir 2 

o)T 


JdK $°(K) JdS r l B^(£) 


(56) 


o . 


Equation (56) can be used to obtain estimates of the residual path deviations 
when the measurements are averaged over many correlation angles. However, 
since wind was also neglected in (56), the results still apply only for 
integration times less than about 1 minute. 



37 


Similar results can be obtained for horizontal propagation paths using 
the parallel path correlation functions. If the ranging site and target are 
stationary, the path separation d is replaced by V,t where is the transverse 
wind velocity. The wind velocity is assumed to be, uniform along the propagation 
path. If the ranging site and target are moving in the same direction at 
identical speeds, is their transverse velocity. The mean square deviation 
of the averaged measurements is 


32 7 r 2 


~/dCB n C€) /dKK$°(K) /dtd-t/T^CKVjt). (57) 


<AL 2 > = 

. avg 

u u u 

When T is large compared to the correlation time L /Vj_ (57) simplifies to 


<AL 2 > - 
avg 


32 it 2 
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(58) 
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3.2 Two Color Ranging Systems 

A significant inT>rovement in ranging accuracy should be possible if two 
laser frequencies are employed instead of one. The difference between the 

A, 

refraction at' different wavelengths can be used to eliminate , to a great extent, 
the effects of atmospheric refraction on the range measurements [13-15]. 

The technique entails measuring the range at two frequencies. The difference 
in range at the two frequencies is a measure of the refractive conditions 
existing over the propagation path at the instant the measurements were taken 
and can be used to improve the accuracy of the range measurements. 

j 

Because the two color system effectively measures conditions along the 
actual propagation path at a given instant-, the system can also correct for 
the random path deviations caused by turbulence. However, dispersion and 
refractivi'ty gradients transverse to the propagation path will cause the two 
color paths to separate slightly. Since the turbulence on the two paths will 
only be partially correlated, path separation reduces the effectiveness of 
the two color system in eliminating errors due to random' fluctuations in 
atmospheric density along the propagation path. This may represent a fundamental 
limitation in the accuracy that can be obtained with two color systems. The 
problem boils down to determining how well the integrated atmospheric density 
along the two color paths is correlated. 

Let L} and L 2 be the measured ranges at wavelengths Aj and A 2 . The actual 
range L is given by [ 13 ] 
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An 12 


AL 12 
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(59) 
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where An i 2 and AL 12 are the differences between the refractive index and 
measured path length. 

Since the two laser beams will traverse slightly different paths in the 
•atmosphere, -the fluctuations on the two paths will be slightly decorrelated. 
Consequently, AL 12 an d AL 2 ^ will have a mean square error of 

<[ALCA 1 )-AL(X 2 )] 2 > = 2 (1-p 12 ) <AL 2 > . (60) 


Pl 2 is the correlation coefficient of the path fluctuations and <AL 2 > is the 
mean square deviation on a single path. 

The correlation coefficient can be calculated from the angular correlation 
function. Young [ 16 ] gives an expression for angular path separation as a 
function of the elevation angle ip and the differential refractivity at the 
ranging site 

\ 

cj) = An 12 cot ip . (61) 


Equation (60) is actually the angular structure function defined in (16) . 

Since the differential path lengths are multiplied by the- coefficient in (S9) , 
the rms error in L caused by partially correlated path fluctuations is given 
by 


L 

error 


n-1 

An 



((j)=An cot \J>) ] 


1 


2 


(62) 


The angular path separation can be decreased by decreasing An, i.e. by 
choosing the two frequencies close together. This decreases the value of the 
structure function and reduces L ' . However, the coefficient (n-1) /An will 

be increased and this tends to increase the error. The dominant behavior of 
L as a function of An can be determined by noting that An cot \p will be 

small compared to the correlation angle . The structure function can be 
approximated by the first term in its power series expansion to give (Sec. 2.3) 
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L - (An)~ l/6 ' (n-1) B 1 '' 2 (“p^)^ (63) 

error <P S 

where 3 is a constant depending on the refractivity spectrum and' C 2 profile. 

The overall error decreases as An increases. It is therefore desirable to 
make An large by choosing a large frequency difference. 

The error reduction provided by a two color system can be estimated from 
equation (62) and the normalized structure functions plotted in Figures 12 and 14 
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error 

AL 
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D (<j>=An cot if/} , 
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(64) 


AL is the rms deviation for a single frequency path. The ratio in equation 
rms 

(64) is plotted in Figure 16 as a function of the elevation angle if). The two 
color system is assumed to use the fundamental and frequency doubled output 
of a laser operating in the near infrared or visible wavelengths -(A j = 1 to .4p, 
*2 = Ai/ 2). Equation (42) (Figure 12) was used for the normalized structure 
function. At the low elevation angles near- 10° the two color system shows 
little improvement over the single color system. For elevation- angles over 
30°, the two color system should provide at least an order, of magnitude 
reduction in the path deviations. 

Similar results would be obtained if the structure functions plotted in 
Figure 14 were used in (64). Figure 14 is based on Brookner's improved C 2 
model (eq. (47)). I'f the Greenwood and Tarazano spectrum is used to calculate 
Dal (curve b in Figure 12) , the ratio in (64) is approximately 1/3 times the 
values given in Figure 16. Thus, the results in Figure 16 probably represent 
a lower bound for the expected improvement of two color systems over single 
color systems. 



Figure 16. Ratio of the rms path deviation for a two-color system 

to the rms deviation for a single color system (eq. (64)) 
ip is the satellite elevation angle. The results were 
calculated for A 1 = 2A„ and apply to the case where 
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For horizontal propagation Thayer [15] calculated the average separation 
between two line-of-sight rays with common endpoints. 
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= Id r A 

12 L W. 


- A 

l dh J 
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(65) 


/ 

where L is the actual range. The error reduction provided by the two color 
system can be estimated using the parallel path structure function 


error 

AL 


rms 


(n-D 

An 


r D AI W=<Ah>) • 

°AlA 


( 66 ) 


Equation (66) is plotted in Figure (17) using the structure function given by 
(34) . The outer scale of turbulence is on the order of 1 to 10 meters for 
propagation near the earth's surface. The two color- system offers significant 
improvement over the shorter paths (<5 km) . However, on the long paths where 
the path deviations can be large, the two color system shows little improve- 
ment over a single color system. 

It should be noted that (66) gives only an approximate value for the 
ratio. The two color paths are not parallel; they are more closely approximated 
by circular. arcs where <Ah> is the average arc separation. In (66) we have 
assumed the arc structure function can be approximated by the parallel path 
structure function evaluated at the average arc separation. 
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L (KM) 

Figure 17. Ratio of the rms path deviation for a two-color system to 
the rms deviation for a single color system (eq. (66)). 

L is the actual range. The results were calculated for 
A^ = 0.6328y and A^ = 0.3164y. 
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4. CONCLUSIONS 

Atmospheric turbulence will introduce errors into distance measurements 
through a number of mechanisms . The most important is the random variations 
in the optical path length. The statistical properties of the path deviations 
were derived in Section 2. The rms path deviations will be typically a few 
millimeters or less. However, when the turbulence is strong (C 2 ~10 -13 m - " 2 / 3 } 
the deviations can be as high as a few centimeters. 

It is possible to reduce the effects of the path deviations by averaging 
many different range measurements. Unfortuantely, the amount of reduction is 
limited by the long tails of the path deviation correlation functions. A 
substantial reduction in the deviations will be obtained only if the measurements 
are averaged over many correlation times (>10, see Sec. 3.1 and Fig. 15}. 

Two color systems can also be used to improve ranging accuracy. The 
difference between refraction at different wavelengths can be used to correct 
the range measurements for the effects of atmospheric refraction and turbulence. 
Over the shorter horizontal propagation paths {<10 km) and at the .higher elevation 
angles (>30°} , the two color systems can reduce the error caused by turbulence 
to a fraction of a millimeter or less. Over the long paths and at low elevation 
angles where the path deviations can be significant, the two color system will not 
perform much better than a single color system. In fact, the two color system 
can have a larger error in some cases (see Fig. 16 and 17} . 

The random path fluctuations caused by turbulence represent a fundamental 
limitation in the accuracy of both single color and two color ranging systems. 

When the turbulence is strong and the effective propagation path long, both systems 
can experience ranging errors approaching a few centimeters. When the turbulence 
is moderate or weak, the errors will be a millimeter or less. 
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APPENDIX A 

Evaluation of D A ^(d) for the Greenwood and Tarazano Spectrum 
The parallel path structure function can be calculated from equation 

(ID 

co 

D AL (d) = a JdK K$°(K) [1-J (Kd)]. (A-l) 

o 

The integral in (A-l) is evaluated by rewritting the Bessel function in 
integral form and substituting the G-T spectrum (eq. (32)) for 4° 

D &L (d) = a' - /d6 JdK K' 5/6 (KL o + 1) “ 11/6 (1 - e lKd COS & ) . fA " 2) 

o o 

The. K integral can be expressed as a confluent hypergeometric function [ 19 ] 
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'2 o o 
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Since J Q (Kd) + o as d “ the normalized structure function is given by 
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When d/L Q is small^U can be expanded in a power series C 20 ] 
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Where 
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APPENDIX B 

Evaluation of the Angular Structure Function 

When the angle <f> is small, the cosine in (16) is approximately one and 
the tangent can be approximated by its argument. In this case the angular 
structure function D AL (<I0 can be written in terms of the parallel path 
structure function D^(d) 
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The power series expansion of B A j j (d)/D A ^( co ) for the G-T spectrum was derived 
in Appendix B 
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The values of the first few coefficients {a m } and (b m ) are calculated in 
Appendix B. After substituting (B-2) into (B-l) we obtain 
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If C^ 2 is described by the exponential model given in (21) , (B-3) becomes 

— = l la F (m + 8/ 3) (<f>/4> ) m + 5/3 + b r(2m + 3) (*/* ) 2(m + ^ ] 

D. t W m=o 111 

AL (8-4) 

L sin ip 

where <j> = . The small 4 asymptote is given by the first term 

s h 

s 

in the series 
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APPENDIX C 

Evaluation of the Angular Correlation Function 
To evaluate the angular correlation function using Brookner's model 
for C 2 (eq. (47)), we start with equation (B-l) in Appendix B and relation 
(33) to obtain 
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The parallel path correlation function calculated using the von "Kantian 
spectrum is given by (25) with y = 11/6 
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For Brookner's model, the integrated contributions of C 2 due to the lower 
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atmosphere and the tropopause are defined as 
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Using equations (C-2) , (C-3) and (47) in (C-l) gives 
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where 4> s and <|> T are the lower atmosphere and tropopause correlation 
angles defined in (49) . The integral in (C-4) can be written in terms of 
a hypergeometric function using the integral relation given in Gradshteyn 
and Ryzik [17] and the transformations given in Abramowitz and Stegun [18]. 
The resulting expression for the angular correlation function is 
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